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On spherically symmetrical accretion in fractal media 
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ABSTRACT 



We use fractional integrals to generalize the description of hydrodynamic accretion in 
fractal media. The fractional continuous medium model allows the generalization of the equa- 
tions of balance of mass density and momentum density. These make it possible to consider 
the general case of spherical hydrodynamic accretion onto a gravitating mass embedded in 
a fractal medium. The general nature of the solution is similar to the "Bondi solution", but 
the accretion rate may vary substantially and the dependence on central mass may change 
significantly depending on dimensionality of the fractal medium. The theory shows consis- 
tency with the observational data and numerical simulation results for the particular case of 
accretion onto pre-main-sequence stars. 
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1 INTRODUCTION 

The interstellar medium (ISM) is believed to have a self- 
similar hierarch ical structure over several orders of magni- 
tude in scale dLarsonl 1 198 ll : iFalgarone. Puget & Peraultl 1 19921 : 
tall" 



iHeithausen et all 19981) . Direct H I absorption observations and 



interstellar scintillation measurements suggest that the structure 
extends down to a scale of 10 AU l Crovisier. Dickey & Kazej 



1 19851 : |L 



anger et 



ail 11994 



Faison et al. 



19981) and possibly even 



to a scale of sub-AU mm et al.l 120051) . However, the latter is 
limited by the spatial resolution of the observations. Hence the 
issue is far from being definite even after observational detec- 
tion of l ower limit of self- similarity scale in some ISM com- 
ponents (Goo dman. Barranco. Wilner & HeverJ fl998l) . Numerous 
theories have attempted to explain the origin, evolution and 
mass distribution of thes e clouds (be gan with the hierarchi- 
cal fragmentation picture, iHovld 19531) and it has been estab- 
lished, from both observations ( Elmegreen & Falgarone 1996) 
and numerical simulations ( Burkert, Bate & Bodenheimer 1997; 



iKlessen. Burkert & Batej [l998; Semel in & Comb es 2000), that the 



interstellar medium has a clumpy hierarc hical self-similar structure 
with a fractal dimension 2.5 < D < 2.7 dSanchez. Alfaro & Perez! 
l2005h in 3 dimensional space. The main reason for this is still not 
properly understood but can result from the underlying fractal ge- 
ometry that may arise due to turbulent processes in the medium. 

Here we have investigated physical processes, in particular hy- 
drodynamics, in such a medium with fractal dimension. We have 
considered the simplest situation of hydrodynamic spherical steady 
accretion of the medium onto a large gravitating mass embedded in 
this medium. We have assumed the dimensionality of the medium 

* E-mail: nirapam@ncra.tifr.res. in 



to be isotropic and homogeneous. It implies that the mass den- 
sity (or, equivalently, dimensionality) is same for any surface inde- 
pendent of the orientation. The famous solution of the continuous 
medium c ase is the "B ondi Solution" for spherical hydrodynamic 
accretion (Bondi 1952). Here we have extended that analysis for a 
medium with fractal dimension D (= 3d) ^ 3. 

To describe physical processes in fract al medium, we have 
used fractional integration and differentiation l lZaslavskvl 2002, and 
references therein). We replace the fractal medium by some contin- 
uous medium and the integrals on the network of frac tal medium 
is approximated by fractional integrals teen et alj2003i) . The inter- 
pretation of fractional integra tion is connected with fractional mass 
dimension (Mandelbrot 1983). Fractional integrals can be consid- 
ered as integrals over fractional dimension space (up to a numerical 
factor; Tarasov 2004). We have chosen the numerical factor prop- 
erly to get the right dimension of any physical parameter and to de- 
rive the standard expression in the limit d — ► 1. The isotropic and 
homogeneous nature of dimensionality is also incorporated prop- 
erly. This model allows us to describe the hydrodynamics in a self- 
consistent way in a fractal medium. 

In this paper we first derive, in §2, the steady state mass accre- 
tion rate for fractal spherical hydrodynamic accretion using sim- 
ple dimensional analysis. In §3 we derive the steady state hydro- 
dynamic equations to describe the spherical accretion in a fractal 
medium and the steady state accretion rate in terms of boundary 
conditions at infinity is derived in §4. We discuss the actual astro- 
physical implication of our analysis and compare our central re- 
sult with observational data and numerical simulation results in §5. 
Possible limitations of our analysis is discussed in §6. Finally, we 
summarize and present our conclusions in §7. 
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2 DIMENSIONAL ANALYSIS OF FRACTAL 
ACCRETION 

We assume a medium that, on a range of length scale R, has a frac- 
tal structure of dimensionality D = 3d < 3 embedded in 3 dimen- 
sional space. Here D refers to the mass dimension and it implies 
that in such a medium of constant density p, the mass enclosed in a 
sphere of radius r will be 



kr 



3(d-l) 



(i) 



where l c is a characteristic inner length of the medium and can take 
arbitrary value in the limit d — > 1. It is the scale below which the 
medium will be continuous. We can define the modified "density" 
for this fractal medium as p = p/l^ d so that Md ~ ~pR 3d . 

For steady state hydrodynamic spherical accretion onto a cen- 
tral mass M from its surroundings fractal medium with a mass di- 
mension of D, the relevant physical parameters will be (1) sound 
speed at a large distance away from the accretor (a^), (2) modified 
"density" of the fractal medium at a large distance away from the 
accretor (p^ = Poo/lc ) and (3) mass of the accretor scaled 
by the gravitational constant (GM). The dimensions of these three 
parameters are 



[a oa ] = [M]°[L] 1 [T] 1 
[- Pao ] = [M] 1 [L]~ D [T]° 
[GM] = [M]°[L] 3 [T]- 2 . 



(2) 
(3) 
(4) 



It is possible to uniquely construct, from these parameters, a quan- 
tity M with dimension [M] 1 [Z/]°[T] -1 . So, from simple dimen- 
sional analysis we get the mass accretion rate 



M 



p^GM) ' 1 ^ 



C- 



GM\d-i 



3(d-l) 



(5) 



Dimension analysis can not be used to fix the dimensionless con- 
stant C in the above equation and does not give a detailed physical 
picture. But we have used the fractional integrals to derived, in a 
more detailed analysis in the following sections, the mass accretion 
rate and found that to be consistent with the accretion rate derived 
from the dimensional analysis. 



3 FRACTIONAL INTEGRALS AND HYDRODYNAMIC 
EQUATIONS 

The integrals o n network of frac tals can be approximated by frac- 
tional integrals dRen et ai1 l2003) and the interpretation of the frac- 
tional integration is connected with fractional mass dimension. We 
consider the fractional integrals as integrals over fractional dimen- 
sion space (up to a numerical factor) and the fractional infinitesimal 
length for a medium with isotropic mass dimension will be given 
by 



dr = -f-rr^dr, d = D/3<l 
It 1 



(6) 



where the constant is chosen to derive the standard expression in the 
limit d 1. Note that the infinitesimal area and volume elements 
in this "fractional continuous" medium of mass dimension D — 3d 
will be 



dA r 



/- 



2(d-l) (p. 



r 2(d-l) 

r sm 



(7) 



dA e 



dA^ = 



dV = 



„2(d-l) 
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1 
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r Hd- 
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,3(d- 
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rdr sin 



rdrdO 



r drsin 



(8) 
(9) 
(10) 



and hence the mass enclosed in a sphere of radius R for constant 
density p will be 

? 3d 

(11) 



, v ;3(d-i) 3d 3 

which will give the standard expression in the limit d — ► 1. 

Below we consider the mass density and momentum density 
balance in such a fractal medium of mass dimension 3d in the case 
of accretion onto a central gravitating mass M for the case of our 
interest (i.e. the steady state hydrodynamic spherical accretion). 

We consider the infinitesimal volume dV in E 3 and the bal- 
ance of mass density or the conservation of mass in that volume in 
the case of steady state spherical accretion will imply 



d [pu(r/l c ) d - 1 dA r ]=0 =► ±{ P ur 3d - 1 ) 



0. 



(12) 



dr dr 
This is the generalized mass density balance equation or the gener- 
alized continuity equation that can be integrated over a surface of 
constant r to give a constant radial accretion rate 

M = ^-pit/^ < ' 1 ~ d V 3d_1 = Constant (independent of r). (13) 
d 1 

In the infinitesimal volume dV in E 3 , the balance of momen- 
tum density will imply 

d 



dt 



(pVdV) = F M + F s 



(14) 



where V is the velocity vector, F M is the total gravitational force 
acting on the mass contained in the infinitesimal volume dV and 
F s is the total surface force due to pressure acting on the surfaces 
bounding the volume dV and are given by 

GMp 



F r M = 



-dV 



(15) 



Ff = e r .(r/lc) d ~ 1 [(pdA r e r ) r — (pdA r e r ) r+dr + 
(pdAgee)e — (pdA e eg)e+do + (pdA^e^)^ — 
(pdA^) (16) 

and, for spherical accretion where p = p(r), it will be reduced to 



F; s = -IfdArdrir/lc 
dr 



(17) 



The total change of radial momentum is given by 

^-ipudV) = i-(pudV) + ^-[pu 2 {r/l c ) d - 1 d~A;]dr. (18) 
dt dt or 

Combining this and equation d 1 2b . in case of steady state we get 

i GMp 



pu — (r/lc) d x dA r dr — — ^-dA r dr{r II. 
dr dr 



-dV(19) 



and this can be simplified to the generalized momentum density 
balance equation 

du 1 dp GM 

u -r + --f + —r = °- (20) 

dr p dr r* 

This equation, using the equation of state p — Kp 1 (1 ^ 7 |) 
and the boundary condition at infinity, can also be integrated to give 



Accretion in fractal media 3 




Figure 1. Velocity profile (u/a-r/ro, ro = GM/a 2 ^) for D = 2.55 and 
7 = 2 a t A = 0.25A C , 1.00A C , 4.00A C 
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(21) 



where a is the sound speed given by a = (dp/dp) 1 / 2 and 0™ is 
the sound speed at infinity. This is a local conservation law and, 
as expected, has exactly the same form as that of the continuous 
medium. 



4 MASS ACCRETION RATE 

We have solved equations H2\ and ( I20t for a smooth, monotonic 
solution without any singularities in the flow. Following the stan- 
dard derivation of the "Bondi solution", we found that there exist 
unique eigen value solution of the problem for a given d and 7 and 
that the solution must pass through a critical point where 

2 



5'(d,7)a 00 , S(d,y) 



(6d- 1) - 3 7 (2d- 1) 



(22) 



The mass accretion rate may not be steady in a small timescale 
because of the "clumpy" structure of the medium but a time average 
mass accretion rate can also be calculated following the standard 
derivation and is given by 



M = 47rA c (d, 7 ) 



3(d-l) 



' V a 2 ) 



= 47rA c (d,7)p 0O a oo (GM/aL) 2 J F , " 3(1 ~ d) (23) 
where F(M, a x , l c ) and A c (d, 7) are dimensionless parameters 
F(M,a x ,l c ) =GM/alol c (24) 

A e (d,7)= (M -^ 1 " M) [^,7)]^- (M - 1) . (25) 

This is consistent with the accretion rate derived from the dimen- 
sional analysis (see Eqn.lO in Sec[2]l. The solution also uniquely 
corresponds to the physically likely situations properly matching 
the boundary conditions (small velocity at large distance and high 
velocity at small r in accretion flow and the opposite in "wind 



flow"). The velocity profile is shown in figure (0 for a particu- 
lar case where D = 3d — 2.55 and 7=2. The general natur e 
of the solution is similar to the "Bondi solution" (Bondi ll952h . 
Note that the smooth, physically selected solution exists only for 
the transonic case with A = A c . It also corresponds to the max- 
imum accretion rate as there exists no meaningful solution either 
for A > A c or for A < A c . 

Note that value of A c is of or der unity. In particular, for d = 1 
and 7 = |, A c = \ ( Bondi 1952) and for 7 = 1, A c is given by 



A c (d, 7 = l) = 



(3d- 



\(X-3d) 



d 2 



(26) 



which is the limit of equation d25t as 7 — > 1, so that A c is continu- 
ous at 7 = 1. 

The most interesting consequence is the dependence of accre- 
tion rate on the factor F(M,a acs ,l c ). In the interstellar medium 
typical ambient sound speed (which goes as s/T) is ~ 10 km s _1 
for a temperature of ~ 10 4 K. In this case the value of F is given 
by 



F(M, Ooo,Ic) = 8.907 (-^)( 



lOkmsN 2 



1 AU \ 
l c ) 



(27) 



This is simply the ratio of the length-scale of "sphere of influence" 
of the accreting object to the scale length of the fractal structure. As 
we are mainly interested in the medium inside the sphere of influ- 
ence, our approach to the problem of accretion in fractal medium 
will be valid only if F > 1. 

In such situations, this factor can change the mass accretion 
rate substantially and it will crucially depend on the fractional di- 
mensionality d. More interestingly, for the cases where this fac- 
tor is important, the accretion rate M will not be proportional 
to M 2 but will be sign ificantly different from that. Fo r example, 
for 2.5 < D < 2.7 dSanchez, Alfaro & Perezl 120051) or equiv- 
alently 0.83 < d < 0.90, the accretion rate M ~ M a where 
1.5 < a < 1.7. 

In the actual astrophysical situation, however, the 
Bondi- Hoyle accretion rate will pr o bably be mo r e rel- 
evant dBondi & HovTel 1 1944 iBondil Il952t lEdgarl |2004 
and references therein). Bondil i 1952 ) proposed the in- 

terpolation formula (corrected u pto a numerical factor by 

Shim a. Matsuda. Takeda & Sawa di Il985h replacing a x by 



(ai 



3 ) 1/2 and Bondi radius r B = GM/a 2 ^ by Bondi-Hoyle 



radius tbh = GM/(a co + v x ) where Vac is the gas velocity 
relative to the star at a large distance. For a medium with density 
and velocity structure, it is not possible to define a unique and 
the present analysis can not be extended to such a situation. But 
it may be extended to a more idealized case where the medium 
has a fractal structure and the accretor is moving in a velocity Voo 
relative to the medium at a large distance. Following a similar 
interpolation method, the accretion rate in this case will be 



M BH = 



47rA e G 2 p 0O M 2 



)3/2 



F 



-3(l-d) 



(28) 



where F ~ tsyi/Ic and A c is a factor of the order of unity. Here 
we have assumed that the gas velocity relative to the star at a large 
distance will modify the accretion rate in fractal medium in a simi- 
lar way that of in the continuous medium. Though a more detailed 
analysis is required to get the Bondi-Hoyle accretion rate and, in 
particular, F in turbulent medium, it is most likely that it will not 
change our main result that the accretion rate will be proportional 
toM 3 ^ 1 . 
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Figure 2. Comparison of data and theory for accretion rates of PMS stars 
and brown dwarfs. The stellar mass is in Mq and the accretion rate is in 

Mqjt- 1 



5 ASTROPHYSICAL IMPLICATIONS 

Our results have important implications for a number of astro- 
physical problems. One of these, for example, is the problem of 
pre-main-sequence (PMS) accretion. Numerical simulation shows 
that approximating the PMS accretion process as Bondi-Hoyle ac- 
cretion leads to ag reement between simulation and observation 
( Padoan et "al] |2005l) . For a solar mass star and for typical sound 
speed flop ~ 0.2 k m s" 1 in ambient molecular cloud filaments 
jPadoan et ai]|2005h . a fractal scale length of ~ 10 AU will be 
smaller than the size of the sphere of influence and hence our ap- 
proach should be valid in this particular case. Here we compare our 
central result with observational data and numerical simulation re- 
sults. For that we have tak en the accretion rates of P MS stars and 
brown dwarfs compiled bylPadoan eT al. (2005) (fromlNatta et al.l 
|2004 IWhite & H illenbran d2004 Muzer olle et alj2004 and refer- 
ences therein) includes all detections but no upper limit. Though 
the the data do not provide any strong support for one theory as op- 
posed to another one due to scatter in the accretion rate attributed 
to (i) an age dependence of the accretion rate, (ii) variation of , 
doo and Voo and (iii) int eraction of accretion flow with jets and out- 
flows on smaller scale I Pad oan etai] l2005). we find that our the- 
ory is consistent with these data. For example, the best fit value of 
ol is 1.49 i 0.13 for log (M Mcr / ¥ B yr ') ^ -10. E ven when 
we include all the compiled data of Padoa net al.l d2005h . the best 
fit value of a is consistent with our theory within 3u error for 
2-5 < D < 2.7. In figure we have shown the data and the best 
fit for log (Maccr/Afp yr" 1 ) ^ -10 Data from lMuzerolle et al.l 
(2004) and White & Hillenbrand 120041) are shown as filled squares 
and empty squares respectively and the rest of the data are shown 
as empty circles. We also found that for higher accretion rate and 
higher central mass, the exponent is significantly different from that 
of the Bondi-Hoyle accretion. For a smaller central mass, the self- 
gravity that we have neglected in our analysis may change the ac- 
cretion rate significantly. 

It is possible to use a high resolution simulation (e.g. 



iKrumholz et~aT]|2006l) to find out the fractional dimensionality and 
the scale length l c by using equation Jilt and counting the num- 
ber of particle n(R) on boxes of different scale length (K). On the 
other hand, comparing the result with existing numerical simula- 
tion results is not very straightforward as in most of the cases it is 
assumed that M accl = AM 2 and the coefficient A is computed 
for the assumed M 2 dependence on accretion rate. But with the 
existing published simulation results, it is still possible to check if 
the variation of this coefficient A within computational uncertainty 
can consistently accommodate our result and we found that the best 
fit of the observational data and our theoretical estim ate are within 
1.5(T scatter around the computed accretion rate of IPadoan et al.l 
(2005). On the other hand, to get a normalized mean accretion rate 
(see IKrumholz et al]|2006l . for details) of 0.1 for a Mach number 
of 5 and D = 2.55, we get, from Eqn.^, r B /l c ~ 17740. For 
a solar mass star and for a typical sound speed ~ 0.2 km s _1 , 
this translates to l c ~ 1.3 AU. Fractal structure at this scale is not 
observationally ruled out. But one should take this as an order of 
magnitude estimate of the scale length as factors like velocity struc- 
ture or magnetic field in real complicated astrophysical situations 
may alter the accretion rate significantly. 

The other potential implication is the black hole accretion. The 
model growth of galactic center black holes assume that the black 
holes accrete at Bondi rate (e.g.. lSpringel. Di Ma tteo & Hernquist 
2005). If the gas around the black hole is fractal in nature, then one 
should rather use the modified accretion rate for fractal medium. 
Here we would like to mention the caveat that there is no defini- 
tive observational evidence against or for the fractal nature of the 
medium in this case. Even if the medium is fractal, the validity of 
our approach, as mentioned in §3, crucially depends on the mass of 
the accretor, the scale length of the fractal and the ambient sound 
speed. 



6 DISCUSSION 

We have derived the steady state hydrodynamic equations to de- 
scribe the spherical accretion in a fractal medium by replacing the 
fractal medium by a "fractional continuous" model. We have de- 
rived the steady state accretion rate in terms of boundary condi- 
tions at infinity in this simplified situation without considering the 
self-gravity of the material for a medium where the dimensionality 
is independent of position and orientation. Magnetic fields, which 
we have not included in these models, will certainly play a major 
role in determining the dynamics. But even without the inclusion 
of magnetic field we have got the following results 

We have found that there exists a unique solution with max- 
imum mass accretion rate and the general natu re of the solution 
is similar to the "Bondi solution" (Bondii l 19521) even for a fractal 
medium with D = 3d < 3. The mass accretion rate, in cases, may 
differ substantially depending on the fractional dimensionality. In 
particular, the accretion rate is likely to be significantly different 
from the "Bondi accretion rate" and will be proportional to M D_1 
in case of accretion in fractional medium with scale length l c very 
different from GM/a 2 ^. 

One limitation of our analysis is that we have not considered 
the self-gravity of the medium. This is justified provided the central 
mass M is very large. In cases where self-gravitation is not negligi- 
ble, it can change the accretion rate significantly. A more important 
limitation is that the present analysis does not account for turbu- 
lent velocity structure. It will certainly play an important role to 
determine the accretion rate but is probably unlikely to change the 
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mass dependence. In a real astrophysical situation with both density 
and velocity structure in the medium, the analytical mass accretion 
rate will hence not be applicable. Neglecting the magnetic fields, as 
mentioned earlier, is the other major limitation of the analysis. The 
presence of magnetic field is more likely to suppress the accretion 
rate and magnetohydrodynamic simulation in fractal medium will 
be required to make a more definitive statement. But the central re- 
sult, that the accretion rate will be proportional to A/ D_1 will not 
be affected by the addition of magnetic fields. 

Our assumption of the equation of state of the form p = Kp 1 
may seems to be a considerable limitation but that is not likely 
to be the case. In ordinary Bondi-Hoyle accretion in continuous 
medium, changing the equation of state changes cha nges the ac- 
cretion rate by a nume rical factor of the order unity dRuffertfl994l ; 
iRuffert & Arnetli 1994) and hence will not change the result much. 
A more detail justification of th e assumption can be found in 
iKrumholz, McKee & Kleij J2006t) . 



7 CONCLUSION 

We have shown that the accretion rate onto a mass embedded in 
a fractal medium may differ, in some cases, significantly from the 
Bondi accretion rate even in the simplest situation. We have used 
the simple model of accretion onto a mass from a fractal medium 
of mass dimensionality D ^ 3 and derived a self-consistent so- 
lution that matches the "Bondi solution" as D — *• 3. The primary 
result of our investigation is that theoretically accretion rate will 
be proportional to M ^ 1 . The observational accretion rate data 
and the numerical simulation for particular astrophysical problem 
of accretion onto PMS stars is consistent with our result. Our find- 
ings suggest that the fractal structure of the medium around the 
accreting mass is playing a major role to determine the accretion 
rate and its dependence on the central mass. The agreement of the 
theoretical prediction with existing numerical simulation implies 
the consistency of the approach. A number of previously published 
numerical and analytical results have not considered this effect and 
may need to be reconsidered. We leave a more detailed treatment of 
the problem, including the effects discussed in this work, the effect 
of self-gravity and stability of fractal accretion, to a future work. 
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